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ABSTRACT 

The results of a next-to-leading order calculation of heavy quark production in longitu- 
dinally polarized photon-nucleon collisions are presented. At cm. energy y/S = 10 GeV, 
for 7 + p — > c + X, cross sections differential in the transverse momentum and rapidity of 
the charmed quark c and the corresponding asymmetries are presented; also, as functions 
of y/S, integrated cross sections, A-factors and the corresponding asymmetries are given. 
Errors in the asymmetries are estimated and the possibility to distinguish between three 
scerarios differing essentially in the polarized gluon distribution is discussed. 
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In spite of the welth of data on polarized deep inelastic scattering, the size and shape of 
the polarized gluon distribution Ag(x) remains a central problem in Spin Physics. Impor- 
tant progress needs experiments on processes with longitudinally polarized initial particles 
dominated by subprocesses with initial gluons. Such a process is 

j+P^Q(Q)+X, (1) 

where Q(Q) denotes heavy quark (antiquark), and, in general, it is dominated by the 
subprocess 

1 + + Q (2) 

Proposals on experiments closely related to, or in, (1) exist in various stages of approval 
[0- 

On the other hand, the importance of determining higher order corrections (HOC) is 
well known. In this Letter we report the essential results of a calculation of the (next-to- 
leading) HOC. 

The subprocesses contributing to (1) follows: 

(A) At leading order (LO, aa s ) §, | 
(Al) The Born subprocess of (2). 

(A2) The resolved 7 via qq — > QQ and gg — > QQ. These involve the polarized pho- 
ton structure functions AF 9 / 7 and AF g / 7 , known only theoretically; hopefully, more 
information will eventually come from experiments on (1). 

(B) At next-to-leading order (NLO, aa 2 s ) 

(Bl) The loop and Bremsstrahlung (Brems, i.e. 7*7 — ► QQg) associated with (2). 
(B2) The subprocess jq — > QQq 

At NLO, scheme independent cross sections can be rigorously obtained only for the sum of 
resolved and direct contributions. 

Notice, that the Abelian part of (Bl) provides the HOC to 

77 -W; (3) 

these HOC have been determined [Q, [5|]. HOC to the process (3) are of interest in them- 
selves: In searches of the Higgs boson (mass uih), for 90 GeV< ran < 2mw, in future 7 — 7 
colliders, the dominant decay mode is H — > bb. With polarized 7' s, the Born contribution 
to the background 77 — > bb is much suppressed; however, due to gluon Brems, HOC have 
an important effect @, ||. 

A calculation of NLO corrections for (1) already exists ||. Our work, however, makes 
use of a different regularization scheme; also, in the treatment of the soft and collinear 
contributions, contrary to which separates them from the hard parts via a cut parameter, 
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we apply more conventional methods [|J. Thus, in view of the importance of (1), we believe 
that an independent calculation is worthwhile. Comparisons, as much as possible, with the 
results of will be also reported. 

Most conveniently, singularities are eliminated by dimensional methods. For polarized 
reactions this requires extension of the Dirac matrix 75 in n = 4 — 2e dimensions. There 
are several schemes for this, and, as in jp, we follow that of dimensional reduction. This 
scheme violates the Ward identity between the vertex and quark self energy functions, but 
care has been taken by introducing a (finite) counterterm, as discussed in 0]. The wave 
function and mass renormalizations are carried on shell 

In the present case charge renormalization is also required. We introduce 



where p, is an arbitrary mass scale entering in n = 4 — 2e dimensions via the change 
g — > gp e . Let Nif = number of light flavors and b = (lliV c — 2Nif) /6. Then we carry 
charge renormalization by introducing the counterterm 

- - (a £ (jjl r ) b - ~A e (m)) (4a) 



e V v ' 3 

where p,R is a regularization mass. In this scheme graphs containing internal loops of the 
heavy quark Q are subtracted out, so that Q is decoupled. This is consistent with parton 
distributions of which the evolution is determined from split functions involving only light 
quarks, as is our case. 

Finally, the gQQ vertex was renormalized via the Slavnov- Taylor identities [jjj]. 

Now, loop contributions to (Bl) were determined via Passarino-Veltman techniques |§; 
and 2 — > 3 parton contributions by going to the cm. (Gottfried- Jackson) frame of Q(Q) 
and final g [0, |J. Also, (B2) was treated by going to the cm. frame of Q(Q) and final 
q. The integrals listed in the Appendices A and C of |§ have been very useful. Some 



remaining integrals are given in an Appendix of [10]. 

Singularities ~ 1/e 2 appearing in the cross sections of (Bl) are cancelled by adding 
loops and Brems. The cancellation of singularities ~ 1/e appearing in (Bl) and (B2) 
required the addition of factorization counterterms. With V\iViiVz the 4-momenta of 7, 
initial parton and observed Q(Q), define 



s = {pi+p2)'\ t = (p 3 -pi) 2 - m 2 , u = (p 3 - p 2 ) 2 - m 2 



Our counterterms are defined in the MS scheme and have the general form: 

. deleter lK(e) „ , . ( m 2 \ £ dan , s 
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where AP a i,{x) split function (in n = 4 dimensions, see below), x proper dimensionless 
variable, M the factorization scale, K(e) a kinematic factor determined from phase-space 
and color and Adas/dt the Born cross section of a 2 — > 2 subprocess with s and either t 
or u replaced by xs and either xt or xu. For ^g — > QQg the 2 — > 2 subprocess is (2); for 
7<f — > <5<5<? one needs two counterterms, one involving (2) and another involving qq — > QQ 
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The resulting finite cross sections are convoluted with polarized parton distributions 



whose evolution is determined by 2- loop anomalous dimensions [12, 13]. On their basis 



several groups have constructed sets of such distributions, differing mainly in the shape 



and size of Ag. We use throughout the NLO sets of one group |L4]]. Also, we use the NLO 
expression of the running coupling constant a s (/j,) with A = 231 MeV and N — Nif + 1 — 4 
flavours. 



The scheme for extending 75 in n dimensions used in |12|, |13| is not dimensional reduc- 
tion, so the addition of certain conversion terms is necessary. The form of the conversion 
terms is easily given in terms of Eq. (5): In n dimensions the split functions have the form: 

API (x, e) = AP ba (x) + eAP^ a (x) ; (6) 

it is AP ba (x) that depends on the scheme. The conversion terms are determined from the 
difference of AP ba (x) in the different schemes. In dimensional reduction AP ba (x) (= P b e a (x)) 
= 0. Refs. [12|, [L3| use the t' Hooft-Veltman scheme, modified so that AP™ q (x,e) = 
P qq (x, s); then AP^ b (x, ) 7^ 0. In terms of these AP^ b (x), the conversion terms are 

dcr conv _ K(0) p£ da B m 

where das/dt the Born cross section of (5) scaled in the same way. 

For the polarized parton distributions we use the sets A, B and C of Ref. 14 and for 
the unpolarized the most recent version CTEQ5 Jl5| . 



Finally, in the absence of any experimental information, as an estimate of the resolved 
7 contributions, we have used the maximal and minimal saturation LO sets of AF 9 / 7 and 
AF 3 / 7 of [|16|], as well as the set of the asymptotic solutions |17|]. In brief, those of fl7| give 



the largest contributions, whereas those of the minimal saturation set give the smallest. 
The analytical calculations were carried with REDUCE and to some extent with FORM. 

Subsequently we present results for Q =c-quark with m = 1.5 GeV only at 

= 10 GeV, relevant to the experiments (a) and (b) of [|TJ . Higher energies and Q =b- 
quark are considered elsewhere |10j. Also, in relation with (4a), we take fiR = fi. 

Fig. 1 presents results related with the differential cross sections Ada/dpT versus xt = 
2pt/VS, where px the transverse momentum of Q; measurement of such cross sections is 
possible in (b) of JT|. Subsequently we denote by Adas/dpr, Ada res /dpT, and Ada iq /dpT 
the contributions to the physical cross section of (Al), (A2), and (B2) correspondingly, and 

by AdajdpT that of the sum (Al), (Bl) and (B2). We use the scale fi = M = (p^ + m 2 ) 1 ^ 2 ; 
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the stability of our results against variations of \i and M is studied in [[UJ. In Fig. 1(a) the 
cross sections Ada /dp T and the corresponding Adas/dpr (denoted by a *) are determined 
for sets A, B and C of |TJj]. The presented Ada res /dpT and Ada iq /dpr correspond to set 
B and the former to the maximal saturation set of Hl6| . In calculating Adas/dpT and 



Ado Tes j dpT we use the NLO sets of 



(a) 



Fig. 1(b) presents the asymmetries 

Ada jdpT 



A L l(pi 



da/dpr 



for sets A, B and C. The resolved 7 contributions have been left out in view of their 
smallness and of the fact that the stage of their present knowledge does not permit a scheme 
independent calculation. Following usual practice in calculating da/dpx we average over 
n — 2 spin degrees of freedom for every incoming boson. Finally, the errors in Fig. 1(b) 
have been estimated from: 

SA * = J^Tm (9) 

In (9) we use unpolarized cross section a integrated over a bin of xt corresponding to 
Ap T = 0.5 GeV and the conditions of 0a] (P B = 80%, P T = 25%, e = 0.014 and L = 2 
fb^ 1 ). Note that the proposal amounts to better conditions and thus smaller 5 All. 
On the basis of our calculations and Fig. 1 we remark the following: 

(i) Defining by Ada ig /dpx the contribution of (Al) and (Bl) we may introduce the 
i^-factors 



K = A^/A d 4^, K, g = A^/A^ (10) 
dpT dpT dpT dpT 



For sets A and B and xt > 0.3, where perturbative QCD is more trustworthy, the 
.ff -factors are > 1 and, in particular K 79 , fairly large. The largeness of K ig is partly 
due to the fact that yS = 10 GeV is fairly low, so a s (n) is fairly large. 

(ii) Ada^qj dpT are found to change little between sets A, B and C; the reason is that the 
valence distributions Am and Ad vary little. Thus Fig. 1(a) presents Ada iq /dpT for 
only one set (B). 

(iii) In general, Ada iq / dpx are significantly smaller than Ada ig / dpx for sets A and B. 
Note that also in unpolarized photoproduction we find that — > QQq contributes 



much less than the NLO 75 — > QQ, in full accord with [18 



(iv) Most importantly, Fig. 1(b) shows that at xt ~ 0.4 one can distinguish sets A and 
C and perhaps also all A, B, C. 



'°)ln this way the effect of the perturbative HOC, as it is reflected e.g. in the magnitude of ^-factors 
(see Eqs. (10)), is made more clear. This is particularly true for polarized reactions, in which the LO and 
NLO Ag differ significantly. 
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Fig. 2 presents results related with the distributions Ada /dY where Y the cm. rapidity 
of Q with respect to the photon. The subsequent notation is equivalent to the above and 
parts (a) and (b) are equivalent to these of Fig. 1. Here we use the scale \i = M = 2m. 
The errors in Fig. 2(b) have been estimated using Eq. (9) with the conditions of ]l|a] and 
the unpolarized cross section integrated over a bin AY = 1. 

On the basis of Fig. 2 we remark: 

(i) At Y > 0, Ada/dY for sets A, B are significantly larger than for C and than 
Ada iq /dY and Ada res /dY. (Fig. 2(a)) 

(ii) Taking into account the errors, Fig. 2(b) suggests that Y m 1.25 ~ 1.5 is the best 
region to distinguish set C from A or B. The region Y < leads to large errors 
because da/dY, like Ada/dY (Fig. 2(a)), is small. 

(iii) Figs 2(a) and 2(b) suggest that integrating the cross sections in the range 1 < Y < 1.5 
offers perhaps the best possibility to distinguish A or B form C. 

Finally, Fig. 3 presents results related with the integrated cross sections Act versus the 
cm. energy J S 1P = VS in a range including y/S = 10 GeV. This quantity will be measured 
in both experiments |l]a] and [06]. Here the scale is again jji — M — 2m. 

To show clearly the NLO effects, we present i^-factors (Fig. 3(a)) in addition to in- 
tegrated cross sections (Fig. 3(b)). At y/S = 10 GeV the ^-factors for all sets exceed 
K — 1. In the range 7 < ^fS < 14 GeV for sets A and B, K are smooth, but for C, K is 
discontinuous due to the vanishing of Aas at ~ 11 GeV. 

Most important is Fig. 3(c), which presents NLO asymmetries. The error at v^S* = 10 
GeV is estimated using in (9) again the conditions of JT|a]. Under these conditions our 
results show that sets A and C can be distinguished, but not sets A and B or B and C. 
Perhaps the proposed SLAC experiment [|TJ>] , which will give results at somewhat lower y/S 
and, as stated, amounts to better conditions, can distinguish also B and C. 



At v5 ~ 10 GeV of importance is the precise knowledge of the charmed quark mass m. 
We understand that recently the uncertainity in m has been somewhat reduced^. Using 
the set B, we find that varying m in the range 1.35 < m < 1.65 changes Act in the range 
27.11 < Act < 13.34 nb and the asymmetry in 6.26% < A LL < 8.88%. 

Ref. [0] works in the scheme of Ref. 19 and presents most of its results using the "stan- 
dard" set of |2(J (GRSVst); for integrated cross sections (Fig. 3), it also presents results 
for sets A and C of ||14||. Also, as stated, they treat their soft and collinear gluon parts by 



the phase space slicing method of |fL8|| , which separates them from the hard gluon parts 
via a cut parameter. First, we have checked that our virtual, soft and collinear contribu- 
tions are in complete agreement (Appendix C of the first reference of 0). Second, with 
respect to numerical comparisons, care is needed in convoluting the LO cross sections with 



( b )We would like to thank A. Despande for this information. 



5 



proper (LO or NLO) parton distributions. Using also GRSVst, and comparing with their 
corrected results, regarding differential cross sections we find differences not exceeding 3%; 
however, regarding asymmetries of integrated cross sections, although in fair agreement, 
our differences are somewhat larger. 

In conclusion, the COMPASS experiment Ha] is expected to provide useful information 
on Ag. This may be said with more emphasis for the proposed SLAC experiment 



We thank I. Bojak for providing us with the results for several quantities we are com- 
paring and for taking his part in doing comparisons. Thanks are also due to G. Bunce, D. 
de Florian, B. Kamal and J. Korner for discussions, to W. Vogelsang for discussions and 
for providing us the sets of |2(J, to P. Bosted for several communications, to A. Despande 
for useful information and remarks and to V. Spanos and G. Veropoulos for participating 
in part of the calculation. 
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(+)This is also the most up to date review of polarized particle reactions 
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FIGURE CAPTIONS 



Fig. 1. Quantities related with the ^-distributions versus x T = 2p T /\/~S: (a) Polarized 
differential cross sections. The LO cross sections are indicated by *. (b) NLO asym- 
metries for sets A, B and C. 

Fig. 2. Quantities related with the rapidity distributions: (a) and (b) as in Fig. 1. 

Fig. 3. Quantities related with the integrated cross sections for sets A, B and C: (a) Factors 
K = A<t/A<7b (b) LO (indicated by *) and NLO cross sections. The presented cross 
section for jq — > QQq corresponds to set B. (c) NLO asymmetries. 
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Fig. 1 



